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We extend our recent work on the two-fluid hydrodynamics of a Bose-condensed gas by including 
collisions involving both condensate and non-condensate atoms. These collisions are essential for 
establishing a state of local thermodynamic equilibrium between the condensate and non-condensate. 
Our theory is more general than the usual Landau two-fluid theory, to which it reduces in the 
appropriate limit, in that it allows one to describe situations in which a state of complete local 
equilibrium between the two components has not been reached. The exchange of atoms between the 
condensate and non-condensate is associated with a new relaxational mode of the gas. 



PACS numbers: 03.75.Fi, 05.30Jp, 67.40. Db 

Recently the authors have given a microscopic deriva- 
tion of the coupled two-fluid hydrodynamic equations for 
a trapped Bose-condensed gas [Q. In the present Letter, 
we report the results of a major extension of the ZGN 
work which takes into account the effects of collisions be- 
tween the condensate and non-condensate atoms. These 
new equations of motion (which we shall call ZGN') al- 
low us to discuss the dynamics when the non-condensate 
atoms are in local thermal equilibrium with each other 
(due to collisions between the excited atoms) but are not 
yet in complete equilibrium with the Bose condensate 
order parameter. This results in the appearance of a 
new relaxational collective mode related to the transfer 
of atoms between the condensate and non-condensate. 
In order to bring out the new physics in a clear fashion, 
these equations are solved for a uniform gas. A more 
complete derivation and discussion is given in Ref. B . 

The non-condensate atoms are described by the distri- 
bution function /(r,p,t), which obeys the kinetic equa- 
tion (we take h = 1 throughout) 



dt m 



VU-V p f = C 12 [f] + C 22 [f} 1 



(1) 



where the effective potential U(r,t) = J7 e xt(r) + 
2g[n c (r,t) + h(r,t)] involves the self-consistent Hartree- 
Fock (HF) mean field. As usual, we treat the interac- 
tions in the s-wave approximation with g = 4ira/m. Here 
n c (r,t) is the condensate density and n(r, t) is the non- 
condensate density given by 



n(r,t) 



dp 



/(r,P,i). 



(2) 



A kinetic equation essentially equivalent to (Q) in the case 
of a uniform system was given in Refs. plP 



The two collision terms in (Q) are given by 
dp 2 f dp 3 



C 22 [f] =4ng< 



(2tt) 3 J (2tt) 3 
-P4)5(e p + e p2 



dp 4 



xS(p + p 2 - p 3 
x [(1 + /)(! + / 2 )/ 3 / 4 - // 2 (1 + / 3 )(1 + h)] 



(3) 



Cl 2 [/] 



4ng 2 n c 



dpi 



dp 2 / dp 3 



(2tt) 3 

x5(mv c + pi - p 2 - P3)S(e c + e Pl - 
x [S(p - Pi) - S(p - p 2 ) - 8{p - p 3 ) 

x[(l + /l)/ 2 /3-/l(l + / 3 )(l + /3)], 



(4) 



with / = f(r,p,t), fi = /(r,p<,t). Eq. (g) takes into 
account the fact that a condensate atom locally has en- 
ergy e c (r, t) = /x c (r, t) + |m»p(r, t) and momentum mv c , 
where the condensate chemical potential /i c and velocity 
v c will be defined in the next paragraph. On the other 
hand, a non-condensate atom locally has the HF energy 
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£p(r, t) = ^— + U(r, t). This particle-like dispersion rela- 
tion means that our analysis is limited to finite temper- 
atures. It follows from this excitation spectrum that in 
the Landau limit (see later) of our microscopic model, the 
condensate density is equal to the superfluid density and 
the non-condensate density is the normal fluid density. 

To complete our microscopic model, we need to have 
an equation of motion for the complex condensate or- 
der parameter Q(r,t) = -y/n c (r, t)e ie ^ r,t \ This equation 
can be rewritten in terms of n c (r, t) and the condensate 
velocity v c = V0(r, f)/m as 



^-+V-(n c v c ) = -r 12 [/] 



m [ — + v c • V ] v c = -V// c , 



(5a) 
(5b) 
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where the condensate chemical potential (in the Thomas- 
Fermi approximation) is given by 



A< c (r, t) = f/extO) + g[n c (r, t) + 2n(r, t)} 



(6) 



The "source" term Ti2[/] in (|5a| ) is defined in terms of 
the C12 collision term in ([|) as 



ri 2 [/] = J ^gC u [/(r )Pl i)]. 



(7) 



We observe that C12 collisions do not conserve the num- 
ber of atoms in the condensate. 

The detailed derivation |2j of Eqs. is based on 

a field-theoretic formulation of an interacting Bose fluid 
with a Bose broken symmetry. Closely related work can 
be found in Refs. Q and ||. The quantum field oper- 
ators V>(r) are s P n t into a condensate ($ = (?/>)) and 
non-condensate part. The key approximation made 
in obtaining (|l|) and (|5|) is the neglect of the anomalous 
pair correlations such as (i/j(r)ip(r)), which we shall refer 
to as the Popov approximation. In M, the C12 collisions 
(involving one condensate atom) were not included and 
thus the source term Ti2 in (5a) was not present. In the 
extended set of ZGN' hydrodynamic equations which fol- 
low from (|l|) and (|J), Ti2 plays a crucial role. Above the 
Bose-Einstein transition temperature (Tbec)j where the 
Bose condensate order parameter vanishes, the kinetic 
equation (|lj) only involves the Uehling-Uhlenbeck colli- 
sion term C22 given by (j^). This kinetic equation has 
been discussed extensively in recent years |?]] . 

In this Letter, we restrict ourselves to the region in 
which C22 collisions are sufficiently rapid to justify the 
assumption that the excited-atom distribution function 
is approximately described by the local equilibrium Bose 
distribution 



/(r,P,t) 



1 



a 0[^:(p-mv n ) 2 +u -fi] _ 1 



(8) 



Here, the temperature parameter (3, normal fluid veloc- 
ity v„, chemical potential fi, and mean field U are all 
functions of r and t. It is important to appreciate that 
the local non-condensate chemical potential p, which ap- 
pears in (H) is distinct from the local condensate chemical 
potential /i c , as defined in (^). One may immediately ver- 
ify that / satisfies C22[/] = 0, independent of the value 
of fi. In contrast, one finds from (^) that, in general, 
Ci2[/] 7^ 0. This means that even if the excited atoms 
are in dynamic local equilibrium described by (|J), the 
source term r^/] in (||) will, in general, be finite. More 
specifically, we have (see also Ref. ||) 



r 12 [/] = -{i 



-^[(i-im(v„-v c ) ! -fi 



T12 



(9) 



where we have introduced a collision time associated with 
the C12 term in (El), 



x<5(mv c + Pi - p 2 - P3)S(s c + s pi - e P2 - e P3 ) . (10) 

We note that Ti2[/] in (j|) vanishes when /2 = /x c + 
^m(v„ — v c ) 2 . However, as we shall see, a state of com- 
plete local equilibrium cannot be treated simply by set- 
ting Ti2 = 0, which was the implicit assumption made 
in deriving the Landau two-fluid equations in earlier 
work |§ . 

With the assumption / ~ /, one can derive hydrody- 
namic equations for the non-condensate by taking mo- 
ments of ([!]) in the usual way. These are the analogue of 
the condensate equations in (^). Linearizing around the 
static thermal equilibrium solutions, these hydrodynamic 
equations are given by 



d5n 
~dt 
_ d6v n 



(11a) 



miiQ- 



&t 



-V5P - SHVUq 

-2gh V(5n + 8n c ) , (lib) 



= - 0^0 V • 6v n - 6v n • VP 
2 

-^gn c oST 12 



(11c) 



One also has 



"I 1 ' ' 1 : ■ : / T^/C'i P. *) = j393/2(z) , (12) 



(2tt) 2 



P(r,t) 



dp P 2 7, ,m 

/» x 3Q A r »P»*) v -0 
(27r) J dm iv„-u 



1 



/3A 3 



(13) 



where z = e^' 1- ' 7 ) is the local fugacity, A = ^Jlit jmk^T 
is the local thermal de Broglie wavelength and g n {z) = 
Y^l 1 z jl n . In static equilibrium (denoted by 0), we of 
course have v„o = v c o = and fi C Q = /to- Thus it follows 
that ri 2 [/°] = 0. The analogous linearized condensate 
equations of motion are 



dSr 



dt 
ddv c 



- V ■ (n c0 (5v c ) - &Ti2, 
-gV{5n c + 2Sh). 



(14a) 
(14b) 



Finally, the source term 5T\2 in these equations can be 
expressed in terms of the fluctuation in the chemical po- 
tential difference /Miff = A ~ Mc ; 



Ponco x 

— "Mdiff , 

r 12 



(15) 



where t® 2 is the equilibrium collision time obtained from 
(pj) with v c = 0, e c — fico and / equal to the absolute 
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equilibrium Bose distribution. We note that adding (11a) 
and (14a) gives the usual continuity equation for the total 
density. 

We now turn to a discussion of our linearized hydrody- 
namic equations given by (|ll^-(|l5|) for a uniform Bose- 
condensed gas (U ext (r) — 0). This means n , n c0 and Pq 
are independent of position. By straightforward calcula- 
tions , one can then reduce our two-fluid equations to 
three coupled equations of motion for the three variables 
5v c , 5v„ and Sfidie- 



m- 



dt 2 



gn c0 V(V ■ 5v< 

i n 



.) + 2gfi V{V ■ 5v n ) 



+2 ff n c0 V(V • 8w c ) - 
d6fj, diS ( 2 



dt 



gn c o 



V-6\ 



V(V • Sv n ) 
2n c0 Pogn c0 

- V • Sv, 



diff 



(16a) 

(16b) 
(16c) 



Here the relaxation time for the chemical potential differ- 
ence (//diff) due to C12 collisions between the condensate 
and non-condensate atoms is given by the expression H 



1 



Pogn-co 



Pp + 2gn n c0 + |-7og»e 



'12 



5 - 
2 



7o Pa 



(17) 



where we have introduced the dimensionless function 
7o = (g(3o/ Ao)gi/2(zo). If we simply omit the terms in- 



volving 5/idiff in ( |!6aD and (16b), we are left with the two 
coupled ZGN equations for 5v n and Sv c given in Ref. ||. 
We see that our new generalized ZGN' equations give 
rise to a coupling between 5v n and 5v c and the local 
variable <5/tdiff, which describes the relative fluctuation 
in the chemical potentials of the two components. This 
is the most important new result in the present Letter. 

It is convenient to solve the three coupled equations in 
( |l6| ) by introducing velocity potentials 5v c = V<^ c and 
5v n = V</>„, and looking for plane wave solutions. We 
obtain from (16c) 



(1 - iur^SfidiS = gncoTfj. I 4>c - ~ <A 



(18) 



Inserting this expression into (16a) and ( |16b| ) gives two 
coupled equations for (j> n and 4> c , 



muj 4> c = gn C Q 1 



+2gh 



mui 4> n = 2gn cQ 



rgj(l - iujT^) 



k (j) n 



(19a) 



Pogn c oT^ n c0 



3^2(1 - iusT^) n 



+ 



5^b „ - 
3n 



9^(1 



"o 



k 2 ^ n . (19b) 



Taking the limit ujt^ 3> 1, one finds that 5/idiff is de- 
coupled from the velocity potentials 4> c ,n and we recover 
the ZGN results g). As expected, in the extreme limit 
ojt^ 1, the effect of C12 collisions are negligible and 
one can simply omit r 12 . 

In the opposite limit ujt^ — > 0, the equations in ( |l9| ) 
yield two phonon-like solutions, W\ 2 = u% 2k, where 
the velocities are given by the roots of the equation 
u A — Au 2 + B = 0. It can be shown that the coefficients 
A and B are in exact agreement with the analogous coef- 
ficients obtained from the usual Landau two-fluid equa- 
tions H9) . The latter theory uses quite different thermo- 
dynamic variables from those used in the present formu- 
lation, and the explicit proof of this equivalence requires 
a lengthy (but straightforward) calculation. It also turns 
out that the first and second sound velocities (tti and 
U2) given by these results (valid for ujt^ <C 1) are nu- 
merically very close to the velocities given by the ZGN 
approximation (valid for lut^ S> 1). The small differences 
involve terms of order g 2 and thus were not picked up in 
the comparison given in Ref. ||. 

The interesting feature of the linearized ZGN' equa- 
tions in ( |l8| ) and ([ljj]) is the existence of a new mode, 
associated with the condensate and non-condensate be- 
ing out of equilibrium (Sfj, c ^ Sft). To a good approx- 
imation (and exact in the k — > limit), it corresponds 
to a mode in which 6v n = (5v c = 0, with a frequency 
given by u> — —i/r^. In the ZGN limit (C12 = 0), this 
reduces to a zero frequency mode. In the Landau limit 
(C12 large), it is a heavily damped relaxational mode. 
According to (|l9|), this equilibration process also gives 
rise to a damping of second sound whose magnitude rel- 
ative to the mode frequency is peaked at ujt^ ~ 1. 

In Fig. 1, we plot the reciprocals of various relaxation 
times involved in our linearized ZGN' equations, as a 
function of temperature. We note that r° 2 goes to zero 
at Tbec and is much smaller than the mean collision time 
expected for a Maxwell-Boltzmann gas. The extra fac- 
tors multiplying l/r° 2 in ( |i~7j ) ensure that l/r M starts to 
decrease as we approach Tbec from below. This is the 
expected "critical slowing down" seen in all second order 
phase transitions involving an order parameter. Our HF 
mean-field approximation is inadequate in the critical re- 
gion close to Tbec 1 and leads to a spurious finite limiting 
value of n c o at Tbec |j] ■ This removes the divergence in 
1/ T i2; which is also why t m in Fig. 1 is finite at Tbec- If 
we simply put n C Q = in our calculations, we would find 
that both r M and l/r^ would diverge at Tbec- 

Our discussion in the present Letter is based on the 
assumption that C22 collisions produce the local equi- 
librium distribution / in (^). As a result, our hydro- 
dynamic equations dTll) do not explicitly depend on a 



3 



relaxation time associated with C^i- This could be in- 
cluded using the standard Chapman-Enskog approach to 
deal with the deviation of / from the local equilibrium 
function / |@. However, an estimate of the collision 
time associated with C22 is given by the scattering out 
term in (0), 



1 



Airg 2 



x<5(pi 



dpi 



dp2 



(2tt) 3 7 (2tt) 3 
P2 ~ P3 ~ P4)^(e Pl 

x/ 1 / 2 1 (l + / 3 1 )(l + /4 )- 



dp s 

(2tt) 3 



dp4 



(20) 



This collision time is plotted in Fig. 1, both above and 
below Tbec- The fact that t® 2 "C t~u at temperatures 
T > 0.8 Tbec is very important, since it allows for the 
possibility that the non-condensate atoms are in local 
equilibrium with each other but not with the conden- 
sate. Above the transition, the value for t^i we obtain 
is in close agreement with the collision time obtained in 
Ref. @. 

These results for the relaxation times are quite inter- 
esting in their own right. The divergence at Tbec (see 
above remarks) is a consequence of the Bosc distribu- 
tion function being used for ff in ([l0|) and (p0|). If a 
Maxwell-Boltzmann distribution function were used, the 
"divergence" shown in Fig. 1 is removed (the importance 
of calculating collision times using the correct Bose dis- 
tribution has also been noted in Ref. ||). Moreover, we 
see that one should not use the classical gas approxima- 
tion for the collision times when determining the cross- 
over from a collisionless (or mean-field) to hydrodynamic 
regimes. The results in Fig. 1 imply that the hydrody- 
namic domain is much easier to reach at finite tempera- 
tures than might have been expected, since the collision 
time can be much smaller than the analogous classical 
gas collision time. 

In summary, starting from a microscopic model, we 
find that the dynamics of a Bose-condensed gas at finite 
temperatures can be divided into three distinct regimes: 

(1) The collisionless regime in which no collision terms 
are included in the kinetic equation (C12 = C22 = 0). 

(2) An intermediate regime in which C22 collisions be- 
tween the excited atoms establish local thermal equilib- 
rium (wr^ < 1) but the C12 collisions do not keep the 
condensate in equilibrium with the non-condensate . 
The relaxation time r M for this equilibration is found to 
be much larger than the collision time t 22 for reaching 
local equilibrium in the non-condensate (see Fig. 1). (3) 
Complete local equilibrium of the condensate and non- 
condensate, which arises in the limit ujt^ <C 1. This is the 
regime conventionally described by the Landau two-fluid 
equations [^]|]]. As stated earlier, the ZGN' equations 
exactly reproduce the results of the Landau equations, 
even though the local dynamical variables used are quite 
different in the two approaches. 

In this Letter, we have only discussed the normal 
modes of our linearized equations in a uniform Bose gas, 



but similar conclusions are obtained for trapped gases [|| . 
Our general equations can also be used to discuss the 
growth and decay of atomic condensates |p|Jl2|, taking 
the dynamics of the non-condensate into account. 
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FIG. 1. Various collision and relaxation times for 
a uniform Bose gas as a function of temperature, for 
gn — 0.1 fee Tbec- The values are normalized to the classi- 
cal collision time rj 1 = naiVoksT /nm) 1 ^ 2 at T = Tbec, 
obtained from (|2(]) using a Maxwell-Boltzmann distribution 
for /. Here, a — 8iva 2 is the atomic cross-section. The tem- 
perature scale changes abruptly above Tbec- The present 
calculations are not valid in the very low temperature region. 



4 



Figure 1 



10 



8 



u 

CO 



=L 4 







"i — i — i — i — i — i — i — i — r 



l/T? 2 ->! 



i/t°4 



l/T 

I 



22 
t 

/', 



I U--f'- I | | | L 



] — i — i — i — i — i — i — i — r 



T 

1/T C , 



J I I I I I I L 



0.2 0.4 0.6 0.8 1 2 3 4 5 6 7 8 9 10 



BEC 



